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ABSTRACT 

We describe a mathematical formalism for determining the mirror shell nodal positions and detector tilts that 
optimize the spatial resolution averaged over a field-of-view for a nested x-ray telescope, assuming known mirror 
segment surface prescriptions and known detector focal surface. The results are expressed in terms of ensemble 
averages over variable combinations of the ray positions and wavevectors in the flat focal plane intersecting the 
optical axis at the nominal on-axis focus, which can be determined by Monte-Carlo ray traces of the individual 
mirror shells. This work is part of our continuing efforts to provide analytical tools to aid in the design process 
for wide-held survey x-ray astronomy missions. 
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1. INTRODUCTION 

Since the pioneering work of Burrows, Burg, and Giaconni, 1 interest within the X-ray astronomy community 
has continued to increase in x-ray survey missions with improved spatial resolution over a wide held-of-view (say 
~30 arcmin). Aspects of a wide-held x-ray imaging survey mission were recently discussed by Murray et al. 2 ' 3 
Extensive reviews of the science and, to a lesser extent, the technology driving such a mission appear in the 
proceedings of the 2009 Bologna WFXT Workshop. 4 There is enormous potential for advances in the studies of 
stars, compact objects, galaxies, active galactic nuclei, galaxy clusters, and cosmology. 

Monte-Carlo simulations 1, 5-9 are currently extensively utilized to determine the design parameters for such 
optics. Since the number of mirror shells per telescope module is typically large (~50 — 100), these procedures 
are computer intensive. While Monte-Carlo simulation will always be part of the design process, this and previous 
papers 10,11 are part of our continuing effort to develop analytic tools to restrict the required search ranges for 
design parameters and to further aid in the design of wide-held x-ray telescopes. 

Here we provide a mathematical formalism for optimizing the placement of mirror shells along the optical axis 
and the tilt of a commonly employed detector configuration at the focus of a x-ray telescope consisting of nested 
mirror shells with known mirror surface prescriptions. We note that the geometric area available is constrained 
and essentially pre-determinecl by the diameter of the launch vehicle faring, the number of telescope modules 
(constrained by the desired FOV and the focal lengths permitted by the launch vehicle faring and available 
extendable optical benches), and the number of mirror shells per module allowed by mass and manufacturing 
constraints. We therefore adopt the spatial resolution averaged over the held-of-view as the figure of merit. 
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2. VARIANCE IN RAY POSITION ON A FOCAL SURFACE 


2.1 Focal surface 

We adopt a coordinate system in which the positive 3-axis points from the detector array to the telescope so 
that on an arbitrary focal surface S 


zs{x,y) = fs(x,y), (1) 

and the optical axis lies along x = y = 0. In the flat plane perpendicular to the optical axis and passing through 

the on-axis focus of mirror shell j, any ray from that shell is completely characterized by its position (xq. y 0 , Zq) 

in that plane, its wavevector there (k x , k y , k z ), and the graze angles with which it reflected from the mirror 
segments of the shell. The coordinates of the ray on S are given by equation (1) and 

Y^J( z S-z 0 ) = x 0 + (jf'j [ fs(x s ,y s ) ~ z 0 ], (2) 

ys = yo + (^j(zs-zo) = yo + [ fs{x s ,y s ) - z 0 ]. (3) 

In practice the shape of the focal surface is highly constrained by detector technology. 

2.2 Single mirror shell 

In terms of the weighted ensemble average values defined in Appendix A.l, the positional variance of the multiply 
reflected rays from a single mirror shell, j, on an arbitrary focal surface S, is 



4s = ^TI [ « 4 >s < Xj >|) + (< y) >s < yj >|) + « z] >s < Zj >|) ]. (4) 

Here n :l is the number of multiply reflected rays from shell j. 

2.3 Nested mirror shells 

In terms of the weighted ensemble average values defined in Appendix A. 2 the variance of the multiply reflected 
rays from a set of J nested shells on S is 

Os = -ff—J [(<x 2 >s - <X >|) + (< y 2 >s - <y >|) + (<z 2 >s - <Z >|) ] (5) 


J 

N = Y, n r 

j = i 

Using the definitions given in Appendices A.l and A. 2, equation (5) can be written in the form 
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Here we have defined 

J J nj 

w = w i = w h fe ’ ( 10 ) 
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where Wj t k is the weight assigned in the ensemble average to the k - th multiply reflected ray from the j-th mirror 
shell (see Appendices A.l and A. 2). When Wj t k = 1, then Wj = nj and W = N. In order to account for 
dependence on energy E, for example for optics with two segments per mirror shell, the natural weight to use 
for doubly reflected rays is the product of the reflectivities from the primary, Rp, and secondary, i?g, mirror 
surfaces: 


Wj,k = Rp{otp,j,k,E) x Rs(ot Sl j,k,E), (11) 

where otpjk and cxs,j,k are the primary and secondary graze angles for the fc-th ray from the j-th mirror pair. 

We see from equations (7) — (9) that <r| is not a simple sum over the variances of the individual shells. We 
rewrite equation (9) for cr| 2 in the forms 


Qij,S = < Xi >s< Xj >S + <Ui >S< Vi >S + < Zi >s< Zj > s 


2 

a S, 2 


2 

a S, 2 


N 


N- 1 
N 

N- 1 


j j 


£ - £ £ (u^) 


Qij,s 


- j=1 

£ m (' 

1=1 


*= 1 1=1 


W 


j j 


\ ' \ ' / a 1 ?: Wj \ 

9ii,s 2^2^ ( u -2 j 

*= 1 1^» 


(12) 

(13) 

(14) 


2.4 Application to an inverted pyramid of detectors 

We now consider a pyramid of four tilted detectors, with apex pointing away from the nested shells, and each 
occupying a quadrant with one corner on the diagonal intersecting the optical axis [see Figure 1)]. In this case, 
the focal surface S = Sp corresponds to the flat, but tilted, surfaces of the four detectors. The tilt angle we 
denote by A Let shell j be displaced along the optical axis so that the apex of the inverted pyramid is a distance 
Szj from the on-axis focus for that shell (in our coordinate system, Szj < 0 if the apex is further from the shell 
midplane than its on-axis focus). Then we have 


a 2 j S = cijfi + 2 bj t o Szj + Cj } o Sz 2 +2 dj t o tant? + 2 eyo Szj tan$ + fjQ tan 2 A (15) 

Here the coefficients ay o, o, Cj, o, dyo, eyo and fj , o are evaluated in the flat plane perpendicular to the optical 
axis and passing through the on-axis focus of shell j. Expressions for each coefficient in the first quadrant in terms 




Figure 1. Illustration of the inverted pyramid detector setup discussed in the text. Each of four detectors lies within a 
numbered quadrant, with gaps between the detector edges. The corner of each detector diagonally opposite from the 
2 -axis is tilted up by an angle The positive 2 -axis points toward the optic and is the optical axis. The point Oj on the 
2 -axis is the on-axis focus for shell j. The apex of the pyramid lies a distance Szj from the point Oj on the optical axis. 
In our coordinate system the case shown satisfies 8zj < 0. Thus the midplane for a two segment optic would be placed a 
distance / — Szj = f + \ Szj | away from the apex, where / is nominal focal length for the shell. 


of the weighted ensemble average for the appropriate combination of ray position vectors and ray wavevectors 
are given in Appendix B.l. 

To second order in the 5z and in tan d, we also have 


Qij,S ~ a ij, 0 + (bijfi SZi + bji o SZj ) + Cjjfi SZj 5Zj 


+ ( djjfi + djifi) tan d + ( e^o Szj + e^o Szj) tan •& 



fjifi ) tan 2 $. 


(16) 


Again the coefficients a,;yo, o> <kj, o> djjfi, ejjfi , and fjjfi are evaluated in the flat plane perpendicular to 
the optical axis and passing through the on-axis focus of shell j. Expressions for each coefficient in the first 



quadrant in terms of the weighted ensemble average for the appropriate combination of ray position vectors and 
ray wavevectors are given in Appendix B.2. 

Using the notation introduced in Appendix B.3, we can express ct Sd compactly as 


cr| D = A 0 + 2 D 0 tanz? + F 0 tan 2 z? + ^ ( 2 Bjfl Szj + Cj t o Szj +2 Ej }0 5zj tanz? ) 

1= i 

(17) 

j j 

EE [ Szi + B jifi Szj) + Cijfi Szi Szj + 5zi + E jit0 Sj) tanz? ] . 
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3. MERIT FUNCTION 

For X-ray survey applications one desires a large effective collecting area over a broad energy range combined 
with good spatial resolution over a wide FOV. The geometric area available is essentially pre-determined by the 
diameter of the launch vehicle faring, the number of desired telescope modules (which are constrained by the 
desired FOV and, in the absence of extendable optical benches, the focal lengths permitted by the launch vehicle 
faring), and the number of mirror shells per module allowed by mass and manufacturing constraints. In our 
work, we have therefore concentrated on optimizing the spatial resolution averaged over the FOV, by minimizing 
the merit function: 


M = (ct|)m 


Sll o d( t> itT 6 d0 w fov(0, 4>) cr|($, (j)) 

4=o ^ fd=n V 6 de w fov(0,4>) 


(18) 


where 9 is the polar off-axis angle for the incident X-rays, (f) is the azimuthal angle for the incident X-rays, 
and wpov(S,4>) is a viewing angle weighting factor. The quantity ag(9,(j)) is the variance in the position of 
multiply reflected rays reaching the focal surface S. We remind the reader that our previous expressions for cr| 
already include weights for the reflectivity of multiply reflected rays from the mirror surfaces of the shell [see 
equation (11)]. 

For the inverted pyramid of detectors described in §2.4, by symmetry the average in Eq. (18) may be restricted 
to (f> £ [0, 7r/4] . This statement neglects any repositioning of the apex of the pyramid to place the on-axis aim 
point on one of detectors. However, such relocation of the on-axis aim point is not necessary when optimizing a 
survey instrument since the best-focus will not be on-axis. 

Substituting equation (17) for the inverted pyramid configuration into equation (18), we find 


{o-% d )m — (Ao)m + 2(U 0 )m tanz? -I- {Fq)m tan 2 zl + ^ ] [ 2(Bj t o) M$Zj + (Cj,o)M<^4 + tanzl ] 
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where ( Q)m denotes an average of the quantity Q over the FOV just as in equation (18). 



4. FINDING THE OPTIMUM VALUES FOR THE OPTIC DISPLACEMENTS AND 

FOR THE DETECTOR TILT 


4.1 Known mirror prescriptions 

When the prescriptions for the mirror surfaces are known in advance, the minimum value for (&s d )m, and hence 
the best values for the optic displacements {Szk, (k = 1, J)} and for the detector tilt angle d, may be found by 
setting the derivatives of < a | >m with respect to each displacement Szk and with respect to tan d equal to 
zero and solving the resulting equations for the Szk and tan A We find 


dp IJm 
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The coefficients (B' k 0 ) M , (E' k 0 ) M and (. E'- k 0 ) M (needed below) are defined in Appendix B.3. These expressions 
are linear in the Sz k and tail'd, and so their simultaneous solution is in principle accessible via well-known 
methods. 

Solving equation (21) for tail'd, we find 


tand = — 


( D 0 )m + ^ {E' k0 ) M &Zk 


fc=i 


/ (Fq)m- 


(22) 


Substituting equation (22) for tail'd into equation (20), and rearranging terms, leads to a system of linear 
equations for the 5zk which we express in linear algebra form: 


(3 ■ 5z = Y, 


(23) 


where the column vectors 5z and Y are given by 


Sz = { Sz k , (k = 1, J) }, (24) 

y = { (Do)m (EIo)m - (Fo)m (B' k> 0 )m, (k = 1, J) } . (25) 

The element of the matrix (3 in the fc-th row and j - th column is given by 

Pkj = S k j [ (F 0 ) M (C k ,o)M - (El-, 0 )m ] ~ (1 ~ faj) [ (■ Fo) M (Cj k ,o)M - ( E' k 0 ) M {E' jkfi ) M ] , (26) 

where 8 k] is the Kronecker delta equal to 1 when k = j and 0 otherwise. Thus we have arrived at a system of 
J linear equations (23) — (26) for the 5z k and an equation (22) for tand in terms of the Sz k . We suggest that 
the system of linear equations (23) — (26) ought to be linearly independent in the sense that the determinant 
| (3 | 0. However, even if this is true in principle, current wide field telescope designs approach 100 closely 

nested mirror shells, so that numerical precision and convergence may be issues for any computer implementation 
of the solution of these equations. 



4.2 Single mirror shell 

For a single mirror shell j with J = j = 1, the conditions (20) and (21) reduce to 

(&i,o)m + tantl (ei i o)M + & z i ( c i,o)m = 0. (27) 


(di,o) M +taxi'& + Szi (e lt0 ) ms = 0, 


(28) 


with solution 


5zi 


(di,o)M(ei,o)M — (&i,o)m(/i,o)m 
(ci,o)m(/i,o)m — (ei,o)M 


tana? 


(6i,o)M(ei,o)M — (ci,o)m(4i,o)m 
(ci,o)m(/i,o)m - ( e i,o)if 
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4.3 Optimizing mirror prescriptions 

When optimization of the mirror prescriptions is desired the above procedure becomes more complex. For 
example consider the case of so-called polynomial x-ray optics 1 assuming the two mirror segment surfaces and 
J mirror shells, for which M — 1 higher order polynomial terms p m ,j,s(z — Zmid)" 1 , with in = (2 ,M), j = (1, J), 
and s = (1, 2), are added to a Wolter I prescription for the mirror segment radius squared. Then derivatives with 
respect to each of the p m ,j,s must be set to zero and simultaneously solved in addition to the conditions for Szj 
and tan 0. Assuming the p m ,j,s are small enough to permit linearization of the new conditions, the linear system 
to be solved now consists of ( M x J + 1 ) equations (including the equation for tan 0) . 


5. CONCLUDING REMARKS 

Monte-Carlo ray traces will continue to play a big role in the design of x-ray telescopes suitable for wide-held 
survey missions. The ultimate goal of our studies is to provide analytical tools to aid in the design process and 
reduce its complexity. In this paper we have provided a mathematical formalism for determining the mirror 
shell nodal positions and detector tilts for a set of nested mirror shells with known surface prescriptions and 
detector focal surface, expressed in terms of ensemble averages over variable combinations of the ray positions 
and wavevectors in the hat focal intersecting the optical axis at the nominal on-axis focus. 
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APPENDIX A. WEIGHTED ENSEMBLE AVERAGES 

A.l Single mirror shell 

For a single mirror shell, designated with the subscript j , we define weighted ensemble average values of the ray 

positions ( x , y, z ) for a set of multiply rehected rays on an arbitrary focal surface S as 


< Xj >s = 

n-j 

- E 

Wj ' 

J k= l 

W s,k %(j,k),S 


< yj >s = 

U 3 

- E 

•' k— 1 
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Here nj is the number of multiply rehected rays from shell j, and {x(j^ s->y(j,k),Si z (j,k),s) are the ( x,y,z ) 
coordinates of the k- th ray on the surface S given by zs = f(xs,ys )• In addition, w St k is the weight assigned 
to the fc-th ray from shell j and 
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We define weighted ensemble average values of the ray wavevectors (k x , k y , k z ), or any combination 
positions and ray wavevectors, in precisely the same way. 

A. 2 Nested mirror shells 

For a set of J nested shells, the corresponding definitions are 
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of ray 


(A. 2-1) 
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Again, we define average values of the ray wavevectors ( k x , k y: k z ) in precisely the same way. 



APPENDIX B. COEFFICIENT DEFINITIONS 

B.l *ls 

We let the subscript 0 denote a quantity evaluated in the flat plane perpendicular to the optical axis and 
intersecting that axis at the on-axis focus. We let the subscript j denote shell j. We let (x,y,z) be the ray 
position vector and (k x , k y , k z ) be the ray wavevector. The coefficients in equation (15) are 

aj , o = ( < x 2 >j, 0 - < x >j fi ) + ( < y 2 > jt0 - < y >f j0 )• (B.l-1) 
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The expression for z(x, y) on the surface of the flat tilted detector varies with quadrant. For the first quadrant, 
we have z = x + y, and therefore in that quadrant 
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On the surface of the detector in the second quadrant we have z = —x + y, in the third quadrant 2 = — x — y, 
and in the fourth quadrant z = x — y. The expressions for dj,o, e 7j o and f 3 ,o are easily modified accordingly. 
However, because we place the apex of the inverted detector pyramid on the optical axis, the configuration is 
symmetric with respect to reflection about the x and y axes and about the straight lines running at 45° with 
respect to those axes. 

We see from equations (B.l-1) — (B.l-11) that to evaluate the coefficients above we need for each shell j the 
values of the weighted ensemble averages 
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B.2 < X* >5< >s + < Vi >s< Vj >S + < Zi >s< zj >s 

To second order in the Sz and in tan$, the coefficients, for the detector in the first quadrant, i: 
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The same comments from the end of Appendix B.l about modifications for the other quadrants and about 
symmetry again apply. We see from equations (B.2-1) — (B.2-9) that to evaluate the coefficients above we need 
for each shell j the values of the weighted ensemble averages 
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Each of these is included in the list given in equation (B.l-12). 
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B.3 Additional notation to provide compact expressions for cr| and its derivatives 

We introduce additional notation in order to reduce the space required for the equations expressing a 2 s and its 
derivatives. We define 
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